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' ^ I Abstract. Much experimental effort is invested these days in fabricating 

pi • nanoelectromechanical systems (NEMS) that are sufficiently small, cold, and clean, 

CJ , so as to approach quantum mechanical behavior as their typical quantum energy scale 

Tiil becomes comparable to that of the ambient thermal energy ksT. Such systems 

will hopefully enable one to observe the quantum behavior of human-made objects, 

^ I and test some of the basic principles of quantum mechanics. Here we expand and 

elaborate on our recent suggestion [Phys. Rev. Lett. 99 (2007) 040404] to exploit 
the nonlinear nature of a nanoresonator in order to observe its transition into the 
CO ' quantum regime. We study this transition for an isolated resonator, as well as one 

[~^ ■ that is coupled to a heat bath at either zero or finite temperature. We argue that by 

^D [ exploiting nonlinearities, quantum dynamics can be probed using technology that is 

almost within reach. Numerical solutions of the equations of motion display the first 
quantum corrections to classical dynamics that appear as the classical-to-quantum 
transition occurs. This provides practical signatures to look for in future experiments 
^ ■ with NEMS resonators. 
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1. Motivation: Mechanical Systems at the Quantum Limit 

This special focus issue of the New Journal of Physics is motivated by the 
growing interest in the study of "Mechanical Systems at the Quantum Limit". 
Nanoelectromechanical systems (NEMS) [H El El HJ [5] offer one of the most natural 
playgrounds for such a study [HI [3, E] • With recent experiments coming within less than 
an order of magnitude from the ability to observe quantum zero-point motion [9|[T0|[TT]. 
ideas about the quantum-to-classical transition (QCT) [121 US [E] may soon become 
experimentally accessible, more than 70 years after Schrodinger described his famous 
cat paradox [T5j. As nanomechanical resonators become smaller, their masses decrease 
and natural frequencies Q increase — exceeding IGHz in recent experiments [161 HTj . 
For such frequencies it is sufficient to cool down to temperatures on the order of 
50mK for the quantum energy hQ to be comparable to the thermal energy fcsT. 
Cooling well below such temperatures should allow one to observe truly quantum 
mechanical phenomena, such as resonances, oscillator number states, superpositions, 
and entanglement [TH [191 [20l [HI [22] , ^^ least for macroscopic objects that are sufficiently 
isolated from their environment. 

Here we expand and elaborate on our recent suggestion [23] to exploit the nonlinear 
nature of NEMS resonators [23] in order to study the transition into the quantum 
domain. Nonlinear behavior of NEMS (and MEMS) resonators is frequently observed 
in experiments [Tl[25l[26l[271[28l[29l[30l[3ll[32l[33], offers interesting theoretical 
challenges [331 ES, [36l [37], and can also be exploited for applications [3HI [39l HOl HI]. 
Here we study a driven nonlinear nanomechanical resonator, also known as the Duffing 
resonator. We choose the system parameters such that the classical Duffing resonator 
is in a regime in which, in the presence of dissipation, it can oscillate at one of two 
different amplitudes, or one of two dynamical steady-states of motion, depending on 
its initial conditions. This dependence on initial conditions has recently been mapped- 
out by Kozinsky et al. ^3] in an experiment with a nanoresonator. Thermally- induced 
switching between the two dynamical steady-states of motion has also been observed in 
a nanoresonator experiment by Aldridge and Cleland |3Qj . 

Naturally, the question arises whether quantum-mechanically the resonator can be 
driven into a superposition of the two possible steady-states of motion, or at least tunnel 
or switch between them at temperatures that are sufficiently low for thermal switching 
to be suppressed, but not necessarily so low to satisfy the constraint fc^T ^ hfl for 
observing full-fledged quantum mechanics. To answer this question we perform two 
separate calculations on the same Duffing resonator, viewing it once as a quantum- 
mechanical object and once as a classical object. With everything else being equal, this 
allows us to contrast the dynamics of a classical resonator with that of its quantum 
clone. We can then search for a regime in which quantum dynamics just begins to 
deviate away from classical dynamics, providing us with practical signatures that can 
be detected experimentally. 

We calculate the dynamics numerically. We start each quantum calculation with 
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a coherent state, which is a minimal wave-packet centered about some point in phase 
space. We start the corresponding classical calculation with an ensemble of initial 
conditions — typically of A^ = 10^ points — drawn from a Gaussian distribution in phase 
space that is identical to the initial quantum-mechanical probability density. We perform 
these calculations for three qualitatively different situations: An isolated resonator with 
no coupling to the environment; a resonator coupled to a heat bath at temperature 
Tenv = 0; and a resonator coupled to a heat bath at a finite temperature T^nv > 0. In all 
cases we display the calculated quantum dynamics in phase space using the quantum 
Wigner function 

W{x,p,t) = — dx'e-^P'^ {x + x'\p{t)\x-x'), (1) 

where p{t) is the usual density operator, and compare it with the time-evolution of 
the classical phase-space density. We remind the reader that the Wigner function is 
not a true probability distribution as it may possess negative values, particularly when 
the quantum state has no classical analog. Nevertheless, it reduces to the quantum 
probability P{x, t)oi observing the system at position x at time t upon integration over 
p, and vice versa]^ 

The paper is organized as follows. Section [2] gives a general discussion of the 
transition from quantum to classical dynamics, and the corresponding limit oi h -^ 0. 
Section [3] describes the classical and quantum mechanical equations of motion for an 
isolated resonator, and section H] discusses the results of our calculation for the isolated 
case. In Section [5] we describe the dynamics of a resonator that is coupled to a heat bath, 
and in Sections [6] and [7] we present the results of our calculation for a resonator, coupled 
to a zero-temperature and a finite-temperature bath respectively. We summarize and 
discuss the significance of our results in Section [HI concluding with some ideas for future 
directions. 

2. The Quantum to Classical Transition (QCT) 

Everyday macroscopic objects behave according to classical physics, expressed in terms 
of Newton's or Hamilton's equations of motion. Since these objects are composed 
of atoms and molecules, one expects this classical behavior to emerge from quantum 
dynamics under a certain limit of high temperature, large masses, and high energies, 
all characterizing the macroscopic world. It is commonly accepted that the transition 
between the quantum and the classical descriptions occurs by letting h go to zero. Of 
course h is actually a non-zero physical constant with the dimensionality of action, so 
whenever one says that "^ goes to zero" one actually means that some dimensionless 
combination of h with a physical quantity characterizing the system — like its ratio to 
the classical action h/S — goes to zero. 

f There exist other possibiUties for describing quantum dynamics in phase space, see for example [42l 
|43l|44], but we find the Wigner function to be among the simplest. 
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One can see in a number of different ways that classical dynamics should be 
obtained from the quantum description by letting h ^ 0. If one expresses the quantum 
mechanical wave function as -ip = ^/pexp[iS/h], then upon letting ^ — *• the Schrodinger 
equation reduces to the classical Hamilton- Jacoby equation for S [l5l Section 6.4]. The 
Feynman path integral approach to quantum mechanics, in which one sums over all 
paths between two points in space and time, reduces to Hamilton's principle of classical 
mechanics, identifying the classical path for which the action satisfies 6S = |46j . 
Finally, one can look at the equation of motion for the quantum Wigner function ([T]). 
This equation can be obtained using the von Neumann equation for the density operator 

Pit), 

P = -^[Hsys, P], (2) 

where Hgyg is the Hamiltonian of the system. For a general potential U{x), which can 
be expanded in a Taylor series, one gets the quantum Liouville equation [47j . 

[dt + ^d. - ^5p) W{x,p,t) = t ^~'l^l^^^^^ (3) 

If one formally sets h = 0, and unless the derivatives of the Wigner function on the 
right-hand side become singular, the right-hand side is equal to zero and one recovers 
the classical Liouville equation for the distribution in phase spacel^l 

It turns out that this naive approach does not work in general, and that one often 
encounters non-analyticities in quantum mechanics as one takes the limit /i — >■ [IE]. For 
an isolated system, the classical Liouville and master equations (see Section 3) violate 
unitarity, as well as the quantum restriction of the density matrix to be semi-positive 
definite. Thus, the limit h —* cannot be connected smoothly to h = 0|§| Therefore, 
as long as the dynamics considered is of an isolated system, the evolution of a quantum 
observable will in general agree with the corresponding classical average only up to a 
finite time, called the Ehrenfest time, tE- For non chaotic systems, it is expected that 
tE will have some power law dependence on h, tE ^ h~ , though the exact value of 6 
may depend on the model [511 [52l |53] . Nevertheless, we observe here that it is difficult 
strictly to define the time scale at which a separation between classical and quantum 
dynamics occurs, mainly because this time scale seems to depend on the observable one 
chooses to measure. Generally, it seems that there is no single time scale, as described 
by Oliveira et al. [HH [55] . 

The contemporary approach to understanding the classical limit of quantum 
mechanics — stemming from early work of Feynman and Vernon [56] and Caldeira and 
Leggett [5^ — is to look at a more realistic model, in which the system of interest is 

I This is also true for a potential containing terms up to quadratic in the position, even for finite Ti. 

Therefore, the quantum and classical dynamics of a harmonic oscillator are identical in the sense of the 

time evolution of phase-space distributions and the expectation values which can be calculated from 

them. 

§ A simple argument is given by Habib et al. j49| to explain this violation. An interesting approach to 

overcome this difficulty is presented, for example, by Gat [50) . 
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coupled in some way to the external world, generically referred to as 'an environment', 
whether it is another system or a measuring device [121 HSl |58]. The interaction with 
the environment must be taken into account, and one distinguishes between two types 
of systems, depending on whether a measurement is performed on the environment or 
not [59l [601 [6T] . If the environment is not observed, then the appropriate description 
of the system is in terms of the reduced density operator, obtained by tracing the 
full density matrix over the variables of the environment. Time evolution is then 
given by a master equation, and the so called weak form of QCT is obtained by 
comparing quantum and classical distributions. This is what we do here, while noting 
in agreement with Habib et al. [52] that the role of the environment is actually two- 
fold. On one hand, it destroys interference patterns with negative values in the Wigner 
function that are incompatible with a classical probability density. On the other 
hand, it causes the classical fine structure to smear over a thermal scale, producing 
a classical distribution without the infinitely-small fine structure, which is incompatible 
with quantum mechanics owing to the finite value of h and the uncertainty principle. 

If the environment is measured, the reduced density matrix depends on the outcome 
of the measurements, and the evolution is said to be conditioned on the observation 
results, as described for example by Habib et al. [61j. This type of evolution can yield 
effectively classical trajectories, and is called the strong form of QCT. We shall analyze 
the latter case, which requires more advanced tools of quantum measurement theory, in 
a future publication. 

3. Isolated System — Method of Calculation 

We consider a nonlinear resonator, such as a doubly-clamped nanomechanical beam [3H 
|32] . nanowire [621 [63], or nanotube [6ll|65], vibrating in its fundamental fiexural mode, 
and thus treated as a single degree of freedom. The resonator is driven, in any of the 
standard NEMS techniques [H [5] , by an external periodic force and is either isolated 
from its environment — a case treated in this section — or coupled to the environment in 
the form of a heat bath at a temperature Tenv — a case that will be treated later. The 
Hamiltonian for such a driven nonlinear resonator, known as the Duffing resonator, is 

1 o "^fi^ _2 ^ ~4 ~ f^ ~ / A\ 

Hsvs = P H X H — X — xF cosut, (4) 

^ 2m^ 2 4 ' ^ ^ 

where m is the effective mass of the resonator, Q its normal frequency, e the strength 

of nonlinearity, F the driving amplitude, and u the driving frequency. Throughout 

the discussion here, and with no loss of generality, we consider the case of a stiffening 

nonlinearity, with e > 0. 

We change from the physical variables, denoted by tildes, to dimensionless variables 

without tildes, by measuring time in units of l/f2, mass in units of m, and length in 

units of JmVt'^ /e. The dimensionless Hamiltonian Hsys is then measured in units of the 

characteristic energy rn^Vf^ /e., 

Ill 
Hsys = -^P^ + ^x^ + jx'^ - xF cos ut, (5) 
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where x, p, and t are now measured in the dimensionless units, and F = F^/e{mVt^)~2 
and uj = Cj/VL are the dimensionless driving amphtude and frequency, respectively. The 
degree to which we approach the quantum domain is indicated by an increasing effective 
value of % as compared with a measure of the classical action S of the system, both stated 
in terms of the scaled units m'^Q^/e. For the case of transverse vibrations of a doubly- 
clamped beam, as shown by Lifshitz and Cross [21], the nonlinear coefficient e is given 
by mVt^ /(P to within a numeral factor of order unity, where d is the width, or diameter, 
of the beam. Thus, the scaled units of action become mQd"^. We shall reinterpret these 
values later in terms of real experimental masses, frequencies, and widths. 

3.1. Classical dynamics - Hamilton and Liouville equations 
Hamilton's equations give the equations of motion for x and p, 

X = p (6) 

p= — X — x"^ + F cosujt. (7) 

For the isolated resonator, the classical Liouville equation for the phase-space 
distribution function f{x,p,t) is 

df 

-Q^ = -{f^Hsys}, (8) 

where {■, ■} denotes the Poisson brackets. Upon substituting the Hamiltonian (jl]) one 
gets (denoting d^ = d/dx etc.) 

dtf = - [pd^ + {-x~x^ + F cos ujt) dp] /, (9) 

where x and p are treated as independent variables. The initial state f{x,p,t = 0) 
is taken to be a Gaussian distribution centered at some [xcPo] and having standard 
deviations of Ax = Jh/2 and Ap = Jh/2, corresponding to an initial quantum coherent 
state with minimum uncertainty AxAp = h/2. 

We can examine the case of a linear resonator simply by omitting the cubic term 
in (I7j). In this case the result is a set of linear ODE's that can be solved analytically. For 
the nonlinear case we integrate the equation of motion (^^ numerically. The solution 
enables us to calculate ensemble averages and make plots of the density in phase space. 
Instead of keeping track of all the trajectories as a function of time, we calculate a 
histogram of the probability density f{x,p,t). We do so by dividing phase space using 
a grid, where each entry of the grid counts the number of trajectories passing through a 
particular square in phase space at any given time. This grid is then used to calculate 
the classical averages. The discrete nature of the grid and the finite number of initial 
conditions (typically 10^) is the reason for the apparent differences between the initial 
classical and quantum states shown in the Figures below. 
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3.2. Quantum dynamics - The Schrodinger equation 

The quantum Hamiltonian operator of the driven Duffing resonator is the same as (jlj), 

Hsys = h(a^a + -] + — ya^ + aj - \ - [a^ + aj F cosut, (10) 

where we have used the ordinary definition of the ladder operators 

[x + ip), (11a) 



a^ = ^i^(x — in). (lift) 

V2h 



We simply solve the corresponding Schrodinger equation numerically by advancing the 
wave function in discrete time steps At, using the interaction-picture time-evolution 
operator. We do so by expanding the wave function in a truncated energy basis \n), 
with n < N. The value of N, which is typically ~ 50, is determined by the requirement 
that the high-energy basis states are never significantly occupied. We use the wave 
function at time t to calculate and then plot the corresponding Wigner function. 

4. Isolated System — Results 

All results in this section are given for a choice of parameters F = 0.01 and u = 1.018, 
which ensures that the resonator is in the bistability region. The results are shown 
in Figure [1] by plotting the Wigner function next to the sampled classical phase-space 
density, for the initial minimal Gaussian wave packet and its later evolution at three 
different times, measured in periods of the drive T = 2ti /uj. We note that the apparent 
discreteness of the classical distributions, as opposed to the smooth quantum Wigner 
functions, stems from the fact that the former are generated from a histogram on a 
discrete grid. The densities are raised to the power 1/4 for better color contrast, where 
blue denotes positive values, and red denotes negative values. A square of area Ti is 
shown at the bottom right corner of the initial state plot to indicate the scale of h. 

We see that for short times the positive (blue) outline of the Wigner function 
resembles that of the classical distribution (Figure [T]d). In addition, a strong interference 
pattern is evident within this outline. In all the snapshots it is evident that the 
Wigner function has a strong positive density where the classical density is also 
large. Nevertheless, two differences between the quantum and classical distributions 
are evident, namely, the strong interference pattern which exists in the Wigner function 
even in the steady state, and the infinitely fine structure which develops in the classical 
distribution. These differences demonstrate that taking the limit of /i — > does not give 
rise to the emergence of classical dynamics from the quantum description, as indicated 
earlier in Section [2l 

We can, however, use the approach proposed by Berry [48j, where the quantum 
probability density is averaged over x and p to account for the limited precision of a 
typical measuring device. If we do the same here, and replace the value of the Wigner 
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(a) Initial coherent state, t = 

0.5 1 1 0.5 1 
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-0.5 



(d) t = 300r 



Figure 1. Isolated Duffing resonator with h — 0.008. Wigner function (left) and 
classical distribution (right) of the initial minimal Gaussian wave packet and its 
evolution at three later times, measured in terms of the drive period T = 2tt/lu. A 
square of area h is shown at the bottom right corner of the initial state plot. The 
functions are scaled by / ^ /^/"^ for better color contrast. 
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0.6 -0.5 



0.5 -0.5 



(a) Original distributions 



(b) Averaged distributions 



Figure 2. Isolated Duffing resonator with h ~ 0.008. (a) Original Wigner function 
and classical distribution (left and right, respectively), shown after 300 periods of the 
drive, as obtained from the calculations shown in Fig. [TJ (b) The same distributions 
averaged using a Gaussian kernel of area 2h. The functions are scaled to the power of 

1/2. 



function W{x,p) with its averaged value over a small area of size 6x6p around the point 
{x,p), the negative and positive parts of the interference pattern cancel each other out. 
If we perform the same averaging also for the classical distribution, the classical fine 
structure is smeared, bringing the two distributions even closer in appearance, as shown 
in Figure [21 

We use the calculated phase-space densities in order to try to estimate the time 
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(b) n = 2,h = 0.032 
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(c) n = l,ri = 0.008 



(d) ?i = 2,n = 0.008 



Figure 3. Relative divergence of quantum and classical expectation values for the 
isolated driven Duffing resonator, calculated using Eq. (|f2p. One cannot conclusively 
estimate the Ehrenfest times Ie for the cases examined here. 



scale tE at which the quantum and classical correspondence is broken. We do so by 
measuring the distance in phase space between the quantum and classical expectation 
values of the first and second moments of x and p as a function of time, scaled by the 
distance of the classical phase-space point from the origin. This measure, 



A" 



\/((^")-^S)' + ((p")-pS)' 



r. 



\2 _^ (rpn\2 



n 



1,2, 



(12) 



gives us the relative divergence of the quantum expectation values from the classical 
ones, where (■) denotes the quantum expectation value, and the subscript d denotes 
the classical ensemble average. The results for A^ and A^ are shown in Fig. [3] on a 
semi-log scale as a function of time, along with their values averaged over one period 
of the drive. For both values of the effective h that are shown it appears difficult to 
determine the Ehrenfest times tE- Even if one insists on making such a determination, 
it seems to depend on the power of x and p being measured, and the resulting estimates 
are inconclusive. Therefore, we cannot determine the Eherenfest times in our model 
with any real certainty. 
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5. Open System — Method of Calculation 

We now couple the Duffing resonator to an external environment by introducing a 
thermal bath. As usual, this coupling introduces dissipation due to loss of energy from 
the resonator to the bath, as well as fluctuations due to the random forces applied by 
the heat bath to the resonator. We use different approaches to add this coupling to 
either the classical or quantum-mechanical resonators, as outlined below. 

5.1. Open classical dynamics - Langevin equation 

In the classical calculation, dissipation is introduced using a Langevin approach, adding 
two terms to the equation of motion for p. The flrst is a velocity dependent friction force 
—'~fx, and the second is a time dependent random force SF{t) acting on the resonator. 
Note that these additional terms are dimensionless and measured using the units that 
were introduced earlier. In particular, 7 is in fact the inverse of the dimensionless 
quality factor Q. The random force is assumed to be a (5-correlated Gaussian white 
noise, satisfying 

{6F{t)) = 0, (13) 

{SF{t)6F{t')) = 2-fkBTenJ{t - t'), (14) 

where ks is Boltzman's constant, thus introducing the bath temperature T^nv The 
equations of motion become 

X = p, (15a) 

p= — X — x^ — 'jp + F coscut + SF(t), (156) 

and are integrated numerically, as in the case of the isolated system. 

5.2. Open quantum dynamics - Master equation 

To describe the evolution of a quantum Duffing resonator that is coupled to an 
environment we can proceed either by using a quantum Langevin approach [12], or by 
coupling the resonator to a bath of harmonic resonators [l2l [66]. The latter approach, 
which we follow here, adds two terms to the Hamiltonian (ITOl) — a Hamiltonian for the 
bath and an interaction Hamiltonian — resulting in 

Htot = Hgys + Hbath + V. (16) 

The von Neumann equation ([2]) for the total density operator ptot, describing the 
resonator and the bath, is 

Ptot = ^[Htot, Ptot], (17) 

where the resonator alone is now described by the reduced density operator, obtained 
by tracing out the bath degrees of freedom, 

Psys{t) = Trbath {Ptot{t)} . (18) 
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We use a standard interaction Hamiltonian [21] of the Caldeira-Leggett [SZ] type in the 
rotating- wave approximation, 

V = J2{KAa^ + K*bla), (19) 

i 

where the Hi are bihnear couphng constants, the bi are annihilation operators acting on 
the bath oscillators, and a is the annihilation operator of the Duffing resonator (1 11 al and 
b). Note that a different choice of coupling other than Eq. (fT9|) may lead to a different 
Master equation. 

By assuming the interaction to be weak, and by employing the Markov 
approximation which assumes that the bath has no memory, we obtain a standard 
master equation [211 1121 EH] , 

1 7 

Psys = -^[Hsys,Psys] " "^(1 + n){a^apsys + Psysa^a - 2apsysa)) 

- -n{aa^psys + Psysdo^ - 2a'^psysa), (20) 

where n = i^e^^i^BTenv _ j^^-i jg ^j-^g Bose-Einstein distribution through which T^nv is 
introduced, 7 = 2-ng{yL)\K(yL)\'^ is interpreted as the damping rate, gi^t) is the density 
of states of the bath, evaluated at the resonator's natural frequency, and k,{VL) is the 
coupling constant at the resonator's natural frequency. We evaluate these two variables 
at the natural frequency fi = 1, because the shift in the resonator's frequency due to 
nonlinearity is assumed small. 

The resulting master equation ( l20l) is of a Lindblad form [42l [58l [67] , which in the 
Schrodinger picture is given by 

p = hH,p] -W{ClCmP + pClCJ+Y.C^pCl, (21) 

lU ^ m m 

where H may differ from Hgys, for example due to a Lamb shift of the resonator's natural 
frequency. The non-unique operators Cm and C^ operate on the Hilbert space of the 
resonator. Identifying the Lindblad form enables us to avoid solving the master equation 
directly — a solution which demands extensive computation. Instead, we use the Monte- 
Carlo wave function method (MCWF) [591168] . which is computationally more efficient. 
It is based on computing many time evolutions of the initial quantum state, so-called 
quantum trajectories, with a non-hermitian Hamiltonian that implements the stochastic 
action of the bath and is derived from Hgyg. Subsequent averaging over the different 
trajectories yields the correct time evolution of the density operator that satisfies the 
master equation (l20l) . We carry out the MCWF numerically using MATLAB. 

6. Open System at T^nv = — Results 

We begin by considering the coupling of the driven Duffing resonator to a zero- 
temperature bath. Again, we choose the values of the different parameters F = 0.015, 
uj = 1.018, and 7 = 0.01, such that the driven nonlinear resonator is operating in its 
bistability regime. For Tenv = 0, in the quantum dynamics the average phonon number 
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n = 0, and in the classical dynamics we expect to see damping without the smearing 
effect of a random force. Thus, the ability of the environment to induce transitions 
between the dynamic states of the resonator is suppressed. We look at a value of 
h = 0.004, as shown in Fig. HI where the initial state is placed within the basin of 
attraction of the large amplitude classical stable solution. The Gaussian classical phase- 
space distribution all flows, as expected, to a single point — the fixed point corresponding 
to the large-amphtude dynamical state. 

In the quantum dynamics, at short times we again see a general positive outline of 
the Wigner function which is similar to the classical density, as for the isolated Duffing 
resonator in Section HI but it quickly deviates from the classical distribution. The 
uncertainly principle prevents it from shrinking to a point as in the classical dynamics. 
More importantly, we clearly observe that the Wigner distribution has substantial weight 
around the state of small-amplitude oscillations, which is inaccessible classically for 
the chosen initial conditions at T^nv = 0. The quantum resonator can still switch 
between the two stable dynamical states, even though the classical resonator cannot. This 
switching takes place either via tunneling — in analogy with the macroscopic quantum 
tunneling between two equilibrium states of a static system [19] — or via quantum 
activation [MUZQ], although at this point it is impossible for us to distinguish between 
these two processes. 




X X 

(a) Initial coherent state, t = 



0.5 -0.5 





(c) t = SOT 



(d) t = 300T 



Figure 4. A driven Duffing resonator coupled to a zero-temperature heat bath with 
scaled U = 0.004 (other parameters are hsted in the text). Wigner functions (left) and 
classical phase-space distributions (right) of the initial minimal Gaussian wave packet — 
located within the basin of attraction of the large-amplitude dynamical state — and its 
evolution at three later times. The large-amplitude fixed point, towards which the 
classical distribution flows, is encircled in Figure (d). A square of area h is shown at 
the bottom right corner of the initial state plot. The functions are scaled by / ^ f^^^ 
for a better color contrast. 
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7. Open System at T^nv > — Results 
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(a) Initial coherent state, t — Q 



0.5 -0.5 
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(c) t = 180T 



(d) t = SOOT 



Figure 5. The driven Duffing resonator coupled to a heat bath at ksTenv — 2?iri 
with Ti — 0.004. Wigner function (left) and classical distribution (right) of an initial 
minimal Gaussian wave packet — straddling the separatrix between the regions in phase 
space flowing to the two stable states — and its evolution at three later times. A square 
of area h is shown at the bottom right corner of the initial state plot. The functions 
are scaled by / ^ f'^^'^ for a better color contrast. 

Figure [5] shows the calculated results for a Duffing resonator with h = 0.004, coupled 
to a heat bath at a finite temperature ksTenv = ^hQ. This temperature is obtained by 
adding a fiuctuating force according to the fiuctuation-dissipation relation (TT4l) . To 
avoid having a random force whose magnitude dominates the dynamics we reduce the 
damping rate to 7 = 0.001, thus requiring a longer time for the resonator to reach its 
final steady state. The remaining parameters are chosen to be F = 0.006 and u = 1.016 
to ensure that the resonator is still operating in the bistability regime. 

Here we choose the initial Gaussian state to straddle the separatrix between regions 
in phase space that flow to the two stable states. The interference pattern that 
develops for short times within the general positive outline in the Wigner function 
is soon destroyed by decoherence, and becomes jittery in space and time. At long 
times, both distributions peak around the two stable states, nevertheless they differ 
significantly. The classical density is tightly locahzed around the two solutions with no 
overlap, indicating that Tenv is too small to induce thermal switching of the classical 
resonator between the two states, as was observed in a recent experiment [3^. The 
Wigner function, on the other hand, is spread out in phase space, indicating that 
h is sufficiently large for the quantum resonator to switch between the two states via 
tunneling or quantum activation [69l [70] . This is demonstrated more clearly in Fig. [6l 
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Figure 6. Probability densities P{x) for the same four states shown in Fig.[5l obtained 
by integrating the phase-space distributions in the p direction. The y axis is scaled 
differently for better visualization. 
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(a) Phase space densities, t — SOOT (b) Probability densities oi x, t = SOOT 

Figure 7. As in Figs, [^d) and [Hd), only that the classical calculation (right) is 
for ksTenv — 17?iJ7, yielding (a) similar phase-space distributions; and (b) similar 
probability densities. Both distributions in (a) are left unsealed. 



which shows the probabihty distributions P{x), that are obtained by integrating the 
phase-space distributions in Fig. [5] in the p direction. 

Further analysis shows that only for temperatures as high as fc^T = IThQ does 
the classical phase-space distribution become as wide as the Wigner function is at 
ksT = 2hQ, as demonstrated in Fig. [71 Thus, in a real experiment, evidence for 
quantum-mechanical dynamics can be demonstrated as long as temperature and other 
sources of noise can be controlled to better than an order of magnitude. Note also 
that minor differences exist between the high-temperature classical distribution and the 
low-temperature quantum distribution {e.g. the relative heights of the two peaks in 
the probability density P{x), as seen in Fig. H^b)), which may assist in distinguishing 
between the two possibilities. 

As we lower the temperature, for example, down to ksTenv = ^^, we see further 
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deviations of the quantum dynamics from the classical one. Figure [H] shows the unsealed 
Wigner functions of the state at t = SOOT for the two temperatures. In both cases the 
function is almost all positive with small negative parts, and are peaked around the 
two classical stable solutions. Nevertheless, one sees that the similarity to the classical 
distribution is evidently higher for the higher temperature. 
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Figure 8. The driven Duffing resonator coupled to a iieat batli at ksTenv = ?i^ and 
at ksTenv = 2?ir2 witli Ti — 0.004. Unsealed Wigner functions of the state at i = SOOT 
are plotted for the two temperatures. As in Fig. [5l the initial Gaussian state straddles 
the separatrix so that the classical distribution splits towards both of the two stable 
dynamical states. With the decrease in temperature one sees further deviation away 
from classical behavior. 



8. Discussion 



We have calculated the classical evolution and the quantum evolution of a driven Duffing 
resonator, whether isolated or coupled to an environment. For the case of an isolated 
resonator — which is less relevant for experiment, yet still interesting theoretically — we 
have confirmed that the quantum and classical evolutions are essentially different, and 
that the transition from one picture to the other does not take place in a simple manner. 
The non-analytic nature of the limit h -^ appears in the form of strong quantum 
interference patterns in the Wigner function, describing the quantum evolution, which 
are never suppressed in the absence of coupling to an environment. In addition, one 
observes an infinitely fine structure in the classical phase-space distribution, which is 
absent from the quantum distribution, owing to the uncertainty principle that smears 
any structure on the scale of h. We have shown that these differences can be resolved 
by introducing a constraint on the experimental resolution. An uncertainty in the 
measurement on the order of h can both smear the fine structure in the classical 
distribution and average out the interference pattern in the quantum Wigner function, 
giving rise to similarly-looking phase-space distributions, as demonstrated in Fig. O 
In all cases investigated here — whether for an isolated resonator or one that is 
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coupled to an environment — we have found that at very early stages of the evolution 
the quantum Wigner function and the classical phase-space distribution agree with each 
other. As can be seen in Figs. [11(b), Hl^b), and[5t^b), at early times the Wigner function 
contains a positive backbone which closely resembles the classical distribution, even 
though a quantum interference pattern develops within this backbone. One expects that 
up to some Ehrenfest time there will be a corresponding agreement between quantum 
observables and classical averages. Nevertheless, even though such agreement is observed 
we could not consistently and conclusively estimate the Ehrenfest time. 

Only after coupling the driven Duffing resonator to an environment at temperature 
Tenv was it possible to search for the regime of interest that we set out to find at the 
beginning of this work. We wanted to find a regime that would enable us to observe the 
first deviations away from classical behavior as we pass through the so-called classical- 
to-quantum transition. We have found this regime by gradually decreasing the physical 
dimensions of the system, or equivalently by increasing the effective value of h, as 
measured by the physical scale mfid"^, set by the parameters of the system. We have 
found an appropriate set of parameters — h ^ 10^'^, ksT^nv = "^fiQ, and Q ~ 10^ with a 
corresponding forcing amplitude that puts the Duffing resonator in its bistable regime — 
for which the quantum dynamics, as it appears in the Wigner function, looks very 
much like classical dynamics, but at the same time shows clear signatures that it is not 
classical. 

We should ask how far we are from reaching such parameters in realistic systems 
with current technology. Reaching quality factors on the order of hundreds to thousands 
is quite routine with standard suspended elastic beam resonators and somewhat 
optimistic with suspended nanotubes, while exciting these resonators into their nonlinear 
regime and observing hysteresis due to bistability is very common. Reaching the 
required effective h ^ 10^^ implies that one needs mVtd'^ k, 10"'^^ m^kg/s. For a NEMS 
beam [3ll [32] with d ~ lOnm and m ^ lO^^^kg, vibrating at fi ~ lO^Hz, we obtain 
mfld'^ ^ 10~^^ m^kg/s, which is 5 orders of magnitude too large. On the other hand, 
with suspended nanowires [621163] or nanotubes [6ll[65] with d ^ Inm and m ^ 10~^^kg, 
also vibrating ai Q ^ lO^Hz, we obtain mfld'^ ~ 10~^^ m^kg/s, which is exactly what 
we want. Note that lOOMHz is quite a high natural frequency for a nanotube. For a 
lower frequency, one would obtain an even better value for the effective h, but this would 
require working at extremely low temperatures, while for Q ^ lOOMHz we would need 
Tenv ~ lOmK. Thus, we believe that with slight improvement over current technology, 
one would be able to see the classical-to-quantum transition described here, by taking 
advantage of the nonlinear nature of doubly-clamped nanoresonators. 

Clearly, the most striking signature of quantum dynamics is the appearance of 
quantum interference with negative regions in the Wigner function. It would be 
wonderful if, at some point in the future, one could directly measure the Wigner function 
and observe these negative regions, as one can already do, for example, in the case of 
trapped atoms [71]. Before such direct measurement is possible, we suggest to perform 
experiments along the lines of the weak from of QCT that we have studied here. In 
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such experiments one should prepare the system in a particular initial point {x,p) in 
phase space, let it evolve, and observe it at some later time t. Repeating the same 
experiment 10,000 times, will allow one to generate a histogram of measured positions 
x{t), similar to the numerical histograms shown in Fig. [6l We emphasize that in this 
kind of experimental protocol, the quantum state of the system — assuming it is indeed 
quantum — is destroyed in the measurement process, a time t after it is initiahzed, but 
otherwise one is free of any considerations of quantum measurement theory. 

It is possible to recognize clear signatures while performing such an experiment, 
implying that a driven Duffing resonator, working in its bistability regime, is behaving 
in accordance with the laws of quantum mechanics. In both cases one should operate in 
a regime of h and T^nv similar to the one used here, in which thermal switching between 
the two dynamical steady-states of the resonator is exponentially suppressed, whereas 
quantum switching is possible. Under such circumstances one can look for the following 
signatures: 

(i) If the initial coherent state is prepared within the basin of attraction of one of 
the two steady-state solutions of the dynamical system, as we did in section [6l the 
classical ensemble should flow as a whole to this solution at long times, while parts 
of the quantum Wigner function should flow to the other steady-state solution. 
This yields a flnite probability of flnding the quantum resonator oscillating in a 
dynamical state that is inaccessible for the classical resonator. 

(ii) Even if the initial coherent state straddles over both basins of attraction, as we 
showed in section [TJ and both distributions split between the two steady-state 
solutions, at long times there is a nonzero probability of flnding the quantum 
resonator in a position between the two steady-state solutions, while classically 
this probability is zero. 

An interesting question that we have not yet addressed here, is whether the quantum 
switching that we see takes place via tunneling or quantum activation, as suggested by 
Marthaler and Dykman [Ml [70]. We will address this question in a future publication 
where we will consider individual quantum trajectories of our system, possibly while 
performing a continuous weak measurement, along the lines of the strong form of QCT. 
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